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Abstract 



The goal of this chapter is to review recent analytical results about the growth of a 
(static) correlation length in glassy systems, and the connection that can be made 
between this length scale and the equilibrium correlation time of its dynamics. The 
definition of such a length scale is first given in a generic setting, including finite- 
dimensional models, along with rigorous bounds linking it to the correlation time. We 
then present some particular cases (finite connectivity mean-field models, and Kac 
limit of finite dimensional systems) where this length can be actually computed. 

0.1 Introduction 

The "Random first order theory" (RFOT) is one of the most widely discussed the- 
ory of glass formation and glassy phenomena in fragile systems. Its origin is based 
upon the observation by Kirkpatrick, Wolynes and Thirumalai (1) that a family of 
abstract long range spin glass models with "one step replica symmetry breaking" 
(1RSB) (2) presents freezing phenomena that in several aspects resemble the observed 
phenomenology in freezing of fragile liquids and other glassy systems. In this per- 
spective long range spin glasses provide a unifying description of glassy phenomena, 
including dynamical aspects (MCT dynamical singularity and aging phenomena) and 
thermodynamical ones (metastability and Kauzmann like entropy crisis), which have 
given rise to predictions verified in simulations and experiments (3; 4). Despite the 
appeal of the resulting picture, it immediately became clear that two problems had to 
be overcome to be able to apply convincingly the theory to supercooled liquids: 1) The 
disordered interactions of spin glasses models is not a realistic microscopic description 
of liquid systems 2) The 1RSB picture is strongly dependent on the long range char- 
acter of the interactions and presents some typical pathologies of mean-field theories. 
Both problems have attracted a lot of attention. 

The first problem has been basically circumvented by the fact that the most sophis- 
ticated Mean-Field theories based on realistic liquid models of particles in interaction 
give back the 1RSB scenario (5). This suggest a high level of universality of glassy 
phenomena that goes from spin models with random interactions to supercooled liq- 
uids, one can thus expect to understand general properties of the latter ones from the 
study of the former ones, which are simpler microscopically. 

The second problem, that is how to effectively take into account the finite range 
character of the interactions, remains in our opinion the main obstacle against an ac- 
complished theory of the glass transition, despite important contributions. The RFOT, 
based on scaling and phenomenological arguments, proposes an intriguing scenario 
about which features of Mean-Field Theory survive in short range systems and what 
of the picture should be modified. However, a theory based on a first principle analysis 
of microscopic models is unfortunately lacking. Understanding the role of finite-range 
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interactions includes both important questions of principle, but not directly related 
to physical observables such as the possibility of an ideal glass transition in systems 
with short range interactions, and questions of direct practical interest as finding a 
theory for activated dynamics and the cross-over between mode coupling behavior and 
activation. 

This chapter reviews some of the theoretical efforts to understand the relation be- 
tween the physics of finite dimensional glassy systems and their mean field description 
and to include finite range interactions in microscopic models of glassy phenomena. A 
crucial point that an accomplished theory of the glass transition should address is the 
growth of correlations that accompanies the increase of relaxation time as tempera- 
ture is decreased. The relaxation time of a supercooled liquid increases dramatically 
upon lowering its temperature, until reaching experimentally accessible timescales at 
the laboratory glass transition temperature T g . In this range of temperature the usual 
static spatial correlations (for instance the structure factor measured in scattering 
experiments) remain essentially the same as the ones of a high temperature liquid. 
These two facts seem contradictory: the physical intuition relates a large correlation 
time to cooperative relaxation mechanisms involving a large volume of the sample 
(which incidentally is an argument in favour of the universality of glassy phenomena, 
the microscopic details being "averaged out" in this case), hence the very strong in- 
crease of the relaxation time around T g should have a trace in spatial correlations. One 
way to address this puzzle is to define dynamical correlation lengths. Convincing ex- 
periments (4), following numerical simulations (6) and dynamical theories (7; 8) have, 
after a long search, demonstrated for the first time growing dynamical correlations (7; 
8; 4). Another way, which we shall review in this contribution, is to define a static 
correlation length through a point-to-set procedure slightly more involved than the 
two-point function underlying the definition of the structure factor. 

A large part of this chapter will be devoted to the discussion of this length and its 
relation with the relaxation time. In Sec. 0.2 we explain in details the definition of the 
point-to-set correlation length and we show how it allows to prove bounds between this 
correlation length and the equilibrium correlation time that agrees with the intuition 
sketched above. We then discuss two classes of disordered spin models of the 1RSB type 
where spatial aspects can be addressed through analytic techniques, namely models on 
diluted random graphs and finite dimensional models in the Kac limit. In the first class 
of models, considered in Sec. 0.3, each spin interacts with a finite number of other spins 
chosen at random. These models can be solved exactly through the cavity method; 
in agreement with the general bounds, at the point of dynamical transition both the 
relaxation time and the correlation length are divergent. The second class, studied in 
Sec. 0.4, consists in genuine finite dimensional models with a tunable interaction range 
ro- In the limit of large r$ one can compute the point-to-set correlation function and 
associated correlation lengths. This leads to a detailed picture of glassy phenomena, 
with a dynamical and static correlation lengths that do not necessarily coincide. 

We hope the style of presentation adopted in this chapter will provide the reader 
with a global view and some mathematical and theoretical tools which should make 
easier the reading of the original, more formal, literature on the subjects we address. 
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Fig. 0.1 Scheme of the thought experiment (9) underlying the definition of the point-to-set 
correlation length. 

0.2 Definition of the point-to-set correlation function and its 
relation to correlation time 

0. 2.1 Heuristic discussion 

In this section we want to introduce the notion of point-to-set correlation functions 
and to show that the correlation length derived from it is relevant for glassy systems, as 
upper and lower bounds on the correlation time can be inferred from this length. The 
presentation will be first done in an informal way, following the thought experiment 
first discussed in (9). We shall then revisit it with more mathematical definitions and 
sketch the results and the methods of proof of (10). 

Let us consider the thought experiment of (9), schematized in Fig. 0.1. In a first 
time one takes a snapshot of an equilibrium configuration of the system under study, 

1. e. for instance the values of the spins, or the positions of particles depending on the 
model investigated. Let us call g_ this first configuration, depicted on the first panel of 
Fig. 0.1. Suppose now that the configuration of the system is frozen to the value it has 
in (7 outside a given volume B around an arbitrary point, for instance the center of the 
system (see middle panel in Fig. 0.1), and that the interior is thermalized in presence 
of this boundary condition. One thus obtains another equilibrium configuration a', 
which is forced to coincide with a outside B. Consider now the following question: how 
similar are a and a' around the center of the system? A precise notion of similarity 
shall be given below, in any case it is natural to expect that the larger the volume B, 
the less similar should g_ and a' be at its center. Indeed the influence of the boundary 
conditions, which force <r' to be very close to g_ when B is small, becomes less and 
less efficient when the boundary is pushed away. This procedure thus allow to define 
a correlation function between a point (the center of the system) and a set of points 
(the boundary of B), hence the name already mentioned. It is understood that in 
the correlation function the similarity measure should be averaged with respect to the 
configurations a and a'. From this function one can further define a correlation length. 
Taking for simplicity B to be a spherical ball of radius £, we shall indeed define the 
correlation length £ c as the minimal radius which brings the point-to-set correlation 
function (i.e. the average measure of similarity of the center of a; and a') below a small 
threshold fixed beforehand. 

The equilibrium correlation time r c of the system can be defined in a similar fashion, 
as the minimal time necessary for the auto-correlation function (average similarity 
measure at the same point, between one equilibrium configuration and the outcome of 
its evolution during a certain amount of time) to drop below a given threshold. It turns 
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out that the intuition discussed in the introduction, namely that large correlation times 
and large correlation lengths are two intertwined phenomena, can be given a precise 
content with these two definitions of t c and i c . Indeed the rigorous proof of (10) 
we shall sketch below implies that £ c < t c < cxp{£J?}, where we have hidden for 
simplicity several constants. The interpretation of these two inequalities might sound 
disappointingly simple. As £ c measures the radius of a correlated region of the system, 
and as for the center of the system to decorrelate it must receive some information 
from the boundary of the correlated region, the lower bound merely states that this 
information cannot propagate faster than ballistically. On the other hand the upper 
bound follows from the fact that the dynamics of the center of the system is weakly 
sensitive to the outside of the correlated zone, hence it should closely resemble the 
dynamics of the ball of radius £ c without its surrounding environment. The latter case 
being the dynamics of a finite system of volume if (where d denotes the dimension of 
the system), its relaxation cannot be slower than exponential in its volume. 

The fact that these two inequalities have natural interpretations does not mean 
they have a trivial content, but rather demonstrates the relevance of the point-to-set 
definition of the correlation length l c . Recall indeed that one of the puzzle of the 
glass phenomenology is the drastic growth of the relaxation time without significant 
traces in the structure factor, hence in the 2-point correlation length. The bounds 
between r c and £ c show that despite this fact the growth of the relaxation time must 
be accompanied by a growth of a (static) correlation length if the latter is appropriately 
defined. 

0.2.2 More precise definition of the correlation function 

We want now to provide the reader with more formal definitions of the quantities 
discussed informally above, before restating with more details the bounds between the 
correlation length and the correlation time. Note that if point-to-sct correlations are 
relatively new in physics they are quite common in the mathematical literature, in par- 
ticular in the context of the so-called tree reconstruction problem, see for instance (12; 
11). 

For the sake of concreteness we shall consider a model of N Ising spins <7j = 
±1, whose global configuration will be denoted g_ = (oi, . . . , ctjv)- F° r a subset S 
of the variable indexes { 1 , . . . , N} we will call g_ s the configuration of the variables 
in S. The energy function (Hamiltonian) is decomposed as E(a) = Ea=i ^afeo)i 
that is a sum of M terms E a , with the a-th term involving a subset denoted da of 
the variables. For instance in a two dimensional square lattice model with nearest 
neighbor interactions there would be an interaction a for each edge of the lattice. 
The definition of E encompasses however more general cases, in particular multi-spin 
interactions involving more than a pair of spins. It can be convenient in such a case to 
represent the network of interactions as a so-called factor graph (13), see Figure 0.2 
for an illustration. Each variable ai, . . . ,apf is associated to a circle vertex, while 
the interactions E\,... ,Em are symbolized with square vertices. An edge is drawn 
between a variable i and an interaction a if and only if E a depends on <7j. On the 
right panel of Fig. 0.2 is drawn a portion of the factor graph corresponding to a square 
lattice model. 
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Fig. 0.2 Left: general example of a factor graph representing an energy 
E(a) = X!a=i EaisLda)- Each circle represents a variable m, each square an interac- 
tion term E a . An edge is drawn between E a and cr; whenever the o-th interaction depends 
on the i-th variable. Right: the case of a square lattice with nearest neighbor interactions. 



We will use in the following the notation d(i,j) for the distance between two 
variables. This will be taken as the graphical distance on the factor graph, that is 
the number of interactions that have to be crossed along a shortest path linking i 
and j. This notion of distance has the virtue of being well-defined for any topology of 
the interaction network, not only for finite-dimensional models. In the latter case the 
graph distance is equivalent to the Euclidean one. 

The Gibbs-Boltzmann probability measure is defined on the space of configurations 
by /i((x) = exp[— /3E(a)]/Z, with the partition function Z ensuring its normalization. 
Angular brackets (•) will be used to denote averages with respect to this law. In order 
to quantify the amount of correlations induced by the Gibbs-Boltzmann probability 
between one variable, say <7i, and a set B of other variables, one can consider a func- 
tion F(g_ B ) which depends only on the values of the variables in B, and compute the 
correlation between the two as {o~iF(a_ B )) — (o~i)(F(<r B )). If for instance <Xj were com- 
pletely independent from the status of the variables in B this quantity would vanish. 
To make contact with the thought experiment explained above one can choose as a 
particular function F the magnetization of the i-th spin conditioned on g_ B , denoted 



This conditional average is defined for an arbitrary function / as 



(/(£'))£ fl =E/(2Wfc 



(0.1) 



where the conditional probability is restricted to configurations which coincide with 
cr R on the variables in B, 




if g_' B = a 
otherwise 



(0.2) 



Z(s.b) ensuring its normalization. With this choice for F we have thus obtained the 
correlation function between a point i and a set B as 
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G(i,B) = (cTiia'^J - (a.,) 2 = X>fe) ^M^Zs)^ - > (0-3) 

a a' 

where the last expression enlightens the connection with the heuristic discussion above: 
the configuration a' corresponds to the second equilibrium configuration drawn con- 
ditioned on the value in B of the first configuration a. 

We would like to emphasize that the configuration er' constructed in the second 
place is as representative as a of the Gibbs-Boltzmann equilibrium measure, in other 
words for any function / of the configuration one has (14) 

((/(2')k B ) = </(£)>■ (0-4) 

This property follows simply from the properties of conditional probability and its 
proof will be omitted. In the following sections we will consider correlations between 
replicas of the system with the same constraint (i.e. g_' and a" are generated indepen- 
dently from fj,(-\a B ) for fixed q_ B and quenched disorder in the energy function E). 
Given any two functions /(<r) and g(cr) one can define the correlation 

C(f,9) = ((f(a')), B ( 9 (a")}, B ) . (0.5) 

Similarly to (0.4) one can show that this coincides with the correlation between ct' and 
a according to 

C f (/ > ff) = (/(2)(fl(£ , )) SB >- (0-6) 

This identity will play an important role in the analysis of the point-to-set function 
in the Kac limit discussed in Sec. 0.4, and we shall call it conditional equilibrium 
condition. 

We come back now to the definition of the point-to-set correlation function, and 
consider more specifically the case where B is the outside of a ball of radius £ around 
i, i.e. the set of variables at distance larger or equal than £ from variable i, B(i,£) = 
{j\d(i,j) > £}■ We shall call G(i,£) = G(i, B(i,£)) the correlation function for this 
geometry. It is intuitively clear that G(i,£) decreases when the radius £ of the ball 
increases, as farther away sites are less correlated with i. One can thus set a small 
threshold e and defines the correlation length for site i as the minimal distance £ 
necessary to make the correlation G(i,£) drops below the threshold s. In formula, 
£ t (e) = min{£\G{i,£) < e}. 

We turn finally to the definition of the correlation time. We shall consider a single 
spin flip dynamics in continuous time, defined through transition rates W(a — > a'). 
The single spin flip assumption means that these rates vanish whenever a and a' differ 
in more than one variable. We assume the rates to verify the detailed balance condition 

n(a)W(a -> a!) = IJ,(a')W(a' -> a) , (0.7) 

which ensures that the Gibbs-Boltzmann probability is stationary under this dynamics. 
Moreover the rate from a configuration a to the configuration where a single variable 
i has been flipped is assumed to depend on a only through the configuration of at and 
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of the variables at unit distance from i . This last condition is obviously fulfilled by the 
usual Monte Carlo dynamics like the Metropolis or heat-bath (Glauber) rules. The 
equilibrium dynamics of the model is defined by these rates and the initial condition 
g_{t = 0) which is drawn from the equilibrium law fi. The average over the initial condi- 
tion and the subsequent evolution shall be denoted again by angular brackets (...). The 
auto-correlation function of variable i is hence defined as Ci{t) = (<7j (0)<7j(t)) — (c^) 2 , 
and we can assign a correlation time tj to this variable as Tj(e) = mm{t\Ci(t) < e}, i.e. 
as the minimal time for the auto-correlation function to drop below a given threshold 
e. To avoid confusion in the following let us emphasize that in the definition of the 
correlation time we consider the dynamics of the whole system; the constraints on the 
finite ball B only appear in the definition of the correlation length. 

0.2.3 Relation between the correlation length and the equilibrium 
correlation time 

Now that the definitions of the correlation time Tj and correlation length t{ have 
been given more precisely we can restate in a more accurate way the bounds derived 
rigorously in (10): 

Cx£i{e')<n{e)< 1 + ffl cp{C , a |B(i,^(e // ))l} . (0-8) 

where \B\ denotes the number of sites in the ball B. In a finite dimensional setting one 
has \B{tj\ ~ £ d , but Eq. (0.8) is valid for any topology of the interaction graph. The 
small thresholds e 1 and e" are functions of e which go to zero when e vanishes, and C1.2 
are numerical constants which depends on the microscopic details of the Hamiltonian 
and of the dynamics. 

Let us make a series of remarks on this result: 

• For simplicity we have explained this relationship between length and time scales 
for a discrete system of Ising spins, with arbitrary interactions. Its extension to 
more general discrete degrees of freedom is simple and was the case considered 
in (10). It is natural to expect that a similar result will hold for particle systems 
evolving in the continuum, for instance by means of a coarse grained occupation 
number field, see e.g. (15). 

• The relationship between li and Tj holds site by site. This is particularly important 
for inhomogeneous systems with quenched disorder, where the correlation lengths 
and times can vary wildly from site to site. 

• In the limit of zero temperature the constant C2 will diverge: the presence of trivial 
"energetic barriers" can lead to very large correlation times without a growing 
correlation length. 

• The lower and upper bound in Eq. (0.8) are widely separated when I grows, 
which could suggest that the bounds are very far from optimal. However it can 
be argued that with the very weak hypotheses made for its derivation, this result 
can only be marginally enhanced (i.e. at the level of the numerical constants and 
with a dynamical exponent z > 2 in the lower bound). Indeed both bounds can 
be saturated in the low/high temperature regimes of some models which enters in 
the range of validity of the result (see below for a discussion of diluted mean-field 
models). 
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• We would like to emphasize the static character of the point-to-set definition of 
the correlation length. Indeed the expression of the correlation function stated 
in Eq. (0.3) involves only equilibrium averages, and does not make any reference 
to the dynamical evolution of the system. Of course if one wants in practice to 
evaluate this function for a system or a model which docs not admit an analyti- 
cal solution, Monte Carlo simulations will probably have to be used to generate 
thcrmalized configurations of the full and constrained systems, yet this dynamics 
is here only a computational tool and not a part of the definition. 

• The dependence of the correlation length and time on the arbitrary threshold 
e may look rather unsatisfactory at first sight. Fortunately this is not a real 
issue for glassy systems: one expects for them a discontinuous behaviour upon 
approaching the glass transition. The spatial and temporal correlation function 
decay to zero in a two-step fashion, with the appearance of a growing plateau. As 
long as e is smaller than the height of the plateau (called Edwards- Anderson or 
non-ergodicity parameter) the asymptotic behaviour of r(e) and i{e) is essentially 
independent on the choice of this threshold. This is however a concern for more 
conventional critical phenomena where the order parameter grows continuously 
at the transition. 

• It should be acknowledged that these bounds do not apply directly for kinetically 
constrained models (KCM) (16; 17). As a matter of fact the equilibrium measure 
of these models is a trivial product measure factorized over the sites, hence the 
correlation length as defined above is always equal to 1. However the dynamical 
rules defining the KCM violates a technical assumption of permissivity which is 
necessary for the derivation of these bounds. 

0.2.4 Determination of the correlation length in finite dimensional 

systems 

In the following subsection we shall give some explanations on how these bounds can 
be proven, then the focus of the rest of the chapter will be put on two families of 
models where analytical computations can be pushed further. Before that we want 
however to mention some papers, mostly numerical, where finite dimensional models 
have been investigated with a perspective somehow related to the point-to-set corre- 
lation function. 

To the best of our knowledge the first papers where glassy systems in confined 
geometry with boundary conditions self-consistently generated according to the pre- 
scription described in the previous sub-sections were in ref. (18). These works investi- 
gated the effect of the presence of a wall or confining geometries on the dynamics of 
Lennard- Jones binary mixtures. 

More recent literature has focused on the thought experiment of (9). In (19) a 
plaquette spin model of a glass was studied analytically and numerically. The effect 
of boundary conditions on finite size subsystems could be analyzed and led to a de- 
termination of a static correlation length along the lines of the thought experiment 
of (9). _ 

A binary mixture of soft-sphere particles (i.e. a fragile glass-former liquid modeled 
microscopically with particles in the continuum) was considered in (15; 21; 20). The 
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Fig. 0.3 Point-to-set correlation function for a binary mixture of soft-sphere particles, 
from (21). When the temperature is lowered the influence of the boundary of a spherical 
cavity persists for larger radius R. 

procedure described above of freezing an equilibrium configuration outside a cavity of 
a given radius and thermalizing its interior was thus implemented with Monte Carlo 
simulations. This led to a demonstration of the growth of the static correlation length 
upon lowering the temperature of the liquid. In Fig. 0.3 we reproduce a result of (21), 
which shows the point-to-set correlation as a function of the radius of the cavity, for 
various temperatures. The reader is also referred to (15; 21; 20) for the interpretation 
of these results along the lines of RFOT and for a critical discussion of the scaling 
exponents of RFOT. 

0.2.5 Tools for the proof 

Couplings. In this subsection we shall sketch in an informal way the method of proof 
of Eq. (0.8), the reader being referred to (10) for the details. We want in particular 
to introduce a very useful probabilistic concept on which the proof of both the upper 
and the lower bounds of (0.8) relies, namely the construction of a coupling between 
stochastic processes (22). Let us first explain what is a coupling on the simplest case 
of two random variables, for instance two biased coins. The first coin X^> takes value 
-Head with probability p±, Tail otherwise, while the probability of Head of the coin 
X^ is p2- A coupling of these two biased coins is a random variable (X^\ X^) which 
can take four values {(H, H), (H, T), (T, H), (T, T)}, such that if one observes only the 
first (resp. second) element of the couple (X( l \X@)), one sees Head occurring with 
probability p\ (resp. p?). Obviously one trivial way to construct such a coupling is 
to take X^ and X^> as independent copies of the original biased coins. The power 
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of the notion of coupling relies in the possibility to introduce a dependency between 
the two elements of the couple, without spoiling the partial (marginal) frequency of 
observations of head in the first or second position. It is indeed easy to realize that 
there exists an infinity of couplings of these two biased coins parametrized by a real 
number. Let us give one explicit example. Suppose without loss of generality that 
Pi > P2, and set the value of the coupling to be 

!{H, H) with probability p 2 , 
(if, T) with probability pi - p 2 , (0.9) 
(T, T) with probability 1 — pi . 

One can easily check that this is indeed a coupling (the marginal probabilities for the 
iJcad of the two coins are respectively p\ and p 2 ), and that, for given values of p\ 
and p2, it minimizes the probability that X^ 1 ' ^ X^ 2 \ For this reason it is termed 
the greedy coupling of the two random variables. In particular if the two coins are 
identical, p\ = p 2 , then X^- 1 ' is always equal to X^ in this coupling. 

Lower bound. This notion of coupling extends naturally to random variables more 
complicated than biased coins, and also to stochastic processes, that is sequences 
of random variables indexed by a time parameter. The proof of the lower bound in 
Eq. (0.8) relies indeed on such a construction that we shall now explain. 

Let us choose a variable i and a positive integer i, and denote as above B (resp. B) 
the outside (resp. inside) of the ball of radius t around i. Consider now the stochastic 
process (£^(*)>£. M)j where two configurations of the same system evolves simul- 
taneously, defined by the following rules: 

• at the initial time t = 0, the two configurations coincide, g_^(0) = = g_, 
with g_ drawn from the equilibrium Gibbs-Boltzmann measure. 

• at any time t > where a variable j attempts an update of its value, both 
configurations arc modified together according to: 

* if j G B, the configuration er*- 2 * 1 (t) is kept unchanged, while the spin o^- (t) is 
flipped or not according to the transition rates W{q_^{t) — s- a'). 

* if j G B. one determines the probability of crj 1 ' '(t) (resp. (J^p {t)) right after 
the update according to W{q_^(t) -t a') (resp. W{a {2 \t) -> a')). Then the 
new values of {a^ (t), (t)) are drawn according to the greedy coupling (as 
defined above) of these two Ising spin random variables. 

Observed separately each element of this coupling corresponds : 

• for o/)(i) to the original dynamics of the whole system. 

• for &} 2 >(t) to the equilibrium dynamics for the inside of the ball B, submitted to a 
time- independent boundary condition in B, with g_ B (t) fixed forever to the value 
a B it has at time t = 0. 

It turns out that the original temporal and point-to-set correlation functions can be 
computed from appropriate averages over this coupling, that we shall still denote (. . .) 
with a slight abuse of notation. Assuming that the equilibrium magnetization (<Tj) 
vanishes to simplify the discussion, one realizes that d(t) = (erf ^O)^ 1 \t)}, while 
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(2) (2) 

G(i, B) is the limit of (af' (O)crf ; (i)) as t goes to infinity. To obtain the lower bound 
of Eq. (0.8) one has to relate in some way the behavior of the spatial and temporal 
correlation functions. From the above observation this translates into a comparison of 
the value of <Ji in the two processes Ss'{t) and u} 2 \t). The idea of the proof is then 
to exploit the properties of the above defined coupling. In fact for "small" times, that 
is of order smaller than £, most probably cr,- 1 ^) = a\ 2 \t). Indeed, at initial times 
the two configurations coincide everywhere in the system. Moreover, when j £ B tries 
to update, if the two configurations coincide on j and its immediate neighborhood, 
then necessarily crj 1 -* = a^p right after the update, because of the use of the greedy 
coupling of the two update probabilities. In other words, disagreement between the 
two configurations can only take birth in the outside B of the ball of radius £, and 
must propagate from the surface of the ball to its center (23; 24). 

Let us rephrase this reasoning with more explicit formulas. We start from a simple 
inequality on the point-to-set correlation function, 

G(i,B) = lim <af (0)af \t)) < (a ? (2) (0)a 4 (2) (t)} , (0.10) 

t— ¥OC 

which holds for any value of t > 0. Indeed the equilibrium auto-correlation functions 
of reversible Markov processes are decreasing in time 1 . To relate the two coupled 
processes let us define the indicator function X(t) = $ CT (i)/ t \ o-( 2 >(t)- We can thus upper- 
bound the spatial correlation function as 

G(i,B) < (^(Ojof >(*)) - (<r?\0)<rW( t )(l X(t))) + (<rf(0)<xf >(t)(l - X(t))) 
<( < 7,[ 1) (0K (1) (t))+2((l-X(t))) 

= Ci(t) + 2V dis (t) , (0.11) 

where Vdis(t) is the probability that the two parts of the coupling disagrees on the 
value of <ii at time t. As explained above disagreement between the two configurations 
of the coupling has to travel from the boundary of B towards i along a a sequence 
of £ adjacent spins on a path from B to i. An upper bound on T'dis(^) can thus be 
obtained by multiplying the number of such paths with the probability that in the 
interval of time [0, t] all £ sites of a given path attempts to update their configurations 
in the right order (from the boundary inwards). The first factor is obviously smaller 
than the maximum connectivity of a spin raised to the power £. Because the times 
where the spin attempts an update form a Poisson process the second factor is the 
probability than a Poisson random variable of average t is greater than £, and this 
last probability is smaller than (et/£) e . Putting these two factors together one obtains 
that "Pdis(i) < (Ct/£) e , where C is a constant which depends on the connectivity of the 
interaction graph. Now if one sets t = Ti(e), the r.h.s. of (0.11) can be made smaller 
than 2e by taking £ larger than some constant multiplied by Tj(e), hence the lower 
bound in Eq. (0.8). 

1 This is a simple property that can be proven using the spectral decomposition of the evolution 
operator W(cr — > a 1 ). 



xvi Abstract 



Upper bound. The upper bound of Eq. (0.8) is obtained by showing that the auto- 
correlation function Ci(t) of the z-th spin is weakly sensitive to the configuration of 
the system out of the ball of radius £i(e). One can thus approximate Cj(i) with the 
value it would have in the system made only of the interior B of the ball of radius 
£i(e). It is natural that the autocorrelation time of the latter cannot grow faster than 
exponentially in its volume. Consider indeed an arbitrary system with n variables, and 
the following coupling between two copies of its dynamics. It is initialized with two 
arbitrary configurations (<l^(0), g_^(0)), and one performs at later times the updates 
of (oj. (t), d^\t)) according to the greedy coupling of the two transition rates. This 
implies that as soon as the two copies coincide, they remain the same for all subsequent 
times. Moreover we assume the dynamics to be "permissive" , that is there is a constant 
K > such that when j attempts an update, irrespectively of the neighborhood of j 
in (cr^, a^ 2 -*), the probability that cry' = cry' after the update is larger than k. This 
condition of permissivity is satisfied at any strictly positive temperature if there are 
no hard constraints in the model (no infinite energy configurations), with possibly k 
vanishing when T — > 0. We can upper-bound the relaxation time by the coalescence 
time of the coupling: if for two arbitrary initial configurations (as different as they 
can be) the coupled dynamics has coalesced at a given time to, hence "forgotten" its 
initial conditions, then the equilibrium dynamics of a single configuration would have 
as well. What remains to be proven is that, with a large probability, this coalescence 
time is not larger than exponential in the number of variables n. This follows from 
considering a particular sequence of update events that brings the two evolving copies 
of the system to coincide, namely that all variables have attempted to update at least 
once, and that their last update brought the variables to the same value in the two 
parts of the coupling. The probability of the last condition is by definition larger than 
K n , hence on time intervals larger than n~ n coalescence is very probable. 

0.3 Computation of the correlation function in mean-field 
(random graph) models 

In this and the next section we shall discuss simplified models, of mean-field nature (in 
a sense that shall be precised) , for which analytical computations of the point-to-set 
correlation function are possible. 

0.3.1 A reminder on mean-field glassy models 

Many researchers agree that important insights on the physics of glassy systems have 
been gained by the study of apparently remote models, namely mean-field spin glasses 
with multi-spin interactions. The paradigmatic example is the so-called p-spin model, 
with p > 3, defined by the Hamiltonian 

E(<t) = - Jii...i,<nr •■■<*, » (0-12) 

l<ii<— <i p <N 

where the J^...,- are independent Gaussian random variables of zero mean and vari- 
ance E[Jj 2 ip ] = 2NP-1 ■ I n tnc Ising version of the model the degrees of freedom arc 
o\ = ±1, while the crj are reals in the spherical case, with the additional constraint 



Computation of the correlation function in mean-field (random graph) models xvii 



of = N. The properties of this model have been extensively studied, we refer the 
reader to the lecture notes (25; 26; 27) and references therein for the references to the 
original works. We shall content ourselves here with a very brief reminder of its salient 
features. From a thermodynamical point of view, the free-energy of this model ex- 
hibits a non-analytic behaviour at the critical temperature T c (also called Kauzmann 
temperature) . However the dynamical properties of the model changes at the dynamic 
transition temperature > T c : the Gibbs measure gets split onto an exponential 
number of pure states, the equilibration time of the system diverges in the thermody- 
namic limit for T < Tj, hence the dynamics exhibit aging (with a non-trivial effective 
temperature) if the system is initialized in a random configuration. The absence of 
any thermodynamic singularity at Td is due to a compensation phenomenon between 
the exponential degeneracy of the pure states (complexity, or configurational entropy) 
and the internal free-energy of the pure states, such that the total free-energy of the 
system in the intermediate phase between and T c coincides with the analytical 
continuation of the high temperature (liquid) phase. In technical words one has a one 
step of replica symmetry breaking (1RSB) behaviour with Parisi breaking parameter 
m = 1. On the contrary at T c the complexity vanishes, and for T < T c the system 
enters a true 1RSB phase with m < 1. 

Of course the main criticism which can be made against the relevance of such a 
model for the study of structural glasses is its completely different microscopic nature, 
in particular because of the quenched disorder present in its definition. The two vari- 
ations of the p-spin model presented in this and in the next section shall not escape 
this criticism, yet they partially cure some of the pathologies of the model and present 
an even richer phenomenology, which motivates the theoretical interest for them. 

0.3.2 Definition of diluted mean-field models 

The model defined by the Hamiltonian (0.12) together with the above stated variance 
of the couplings constants J^...^ will be termed the fully-connected p-spin model. 
Indeed all of the (^) possible interactions between p-uplets of spins are present in the 
system, even if each of them is individually weak in order to have an extensive scaling 
of the Hamiltonian in the thermodynamic limit. In this case each spin interacts with 
all the others, in other words the distance d(i,j) between any pair of spins defined in 
Sec. 0.2.2 is equal to one. This peculiarity, very different from the finite-dimensional 
models with short-range interactions where any degree of freedom interacts only with 
a finite number of neighbors, can be cured with the following definition of the coupling 
constants: 

j = J±l with probability a ^ 

1 with probability 1 — aj^ 

This case will be called the diluted p-spin model. At variance with the fully-connected 
case, only an extensive number (concentrated around its average aN) of interactions 
are now present, but each of them is strong (of order 1). Moreover a short computation 
reveals that the degree of one spin, that is the number of interactions it belongs to, 
converges in the thermodynamic limit to a Poisson random variable of average op, 
hence finite with respect to the system size N. 



xviii Abstract 



T 




(V 



Fig. 0.4 Shape of the phase diagram for the diluted p-spin model. 



This kind of diluted model, pioneered in the case p = 2 by Viana and Bray (28), can 
be viewed as intermediary between fully-connected and finite-dimensional ones. They 
share with the latter the finiteness of the connectivity of each degree of freedom, but are 
mean-field as the former ones. Indeed there is no a priori underlying Euclidean space in 
their definition. This allows for detailed analytical computation of their properties, and 
provides a family of models where the Bethe-Peierls approximation is actually exact 
in the thermodynamic limit. Moreover there is one more control parameter besides 
the temperature, namely the "density" parameter a which controls the number of 
interactions present in the system. 

There has been an important research effort in the last decades in the community of 
statistical mechanics of disordered systems to develop analytical tools for tackling such 
diluted spin-glasses models. An important motivation has been the intimate connection 
between these models and random combinatorial optimization problems, as the k- 
satisfiability and the coloring of random graphs. An extensive account of this line of 
research can be found in the recent book (29). 

Before turning to some explanations on the computations in diluted models, in 
particular the determination of the point-to-set correlation function, let us plot in 
Fig. 0.4 the shape of the phase diagram in the (density a, temperature T) plane. The 
two transition temperatures and T c of the fully-connected model (which corresponds 
to the a — > oo limit up to a rescaling of the energy) becomes two lines of transitions, 
which ends up at zero temperature for finite values of the connectivity a, respectively 
ay and a c . The zero-temperature properties of this model have been largely studied in 
the context of optimization problems (where it is known as XORSAT), the thresholds 
ay and a c corresponding there to the clustering and satisfiability transitions (30; 31). 

0.3.3 Computation of the correlation function 

We shall now explain the main steps in the computation of the point-to-set correlation 
function in mean-field diluted systems. Let us first discuss the properties of the (ran- 
dom) factor graph associated with the diluted p-spin model defined by Eqs. (0.12,0.13). 
To make contact with the general definitions we rewrite it as 





(0.14) 



a 
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Fig. 0.5 Pictorial representation of the recursive computation on a tree. Left: the effective 
magnetic field Hi is obtained from the effective fields in the sub-trees F a -fi, see Eq. (0.16). 
Right: the construction of the field for Fi_> from the as shown in the second part of 

Eq. (0.17). 

where we have only retained the p-uplets with a non-vanishing coupling constant and 
renumbered them in terms of the interaction index a. Each interaction a involves p 
variables, while each variable appears in a random number of interactions, which is 
easily found to be a Poisson random variable of average ap. What is slightly less 
obvious is the fact that these random factor graphs have a local tree structure, that 
is short loops are relatively rare, and only loops of length of order logiV begin to 
proliferate in the thermodynamic limit (N — > oo with a finite). It is this locally tree- 
like property which allows for analytical computations in these models; note that this 
property is shared by other definitions of the random factor graph model, in particular 
with degree distributions other than Poissonian, for instance deterministic. 

Let us forget for a moment the existence of (long) loops in the factor graph and 
see how the problem can be solved easily on a finite tree. The key point is the natural 
recursive structure of a tree: removing a vertex from a tree leads to a set of smaller sub- 
trees, which, by definition, arc disconnected one from the others. So if one can solve 
the considered problem on trees of some given size, and glue together the solutions to 
compute the solution on the concatenated tree, then by recurrence the solution can be 
found for any tree. Consider for instance the computation of the local magnetization 
(ex;) of a spin in a tree factor graph. According to the recursive strategy just explained, 
this can be expressed in terms of the magnetization in the sub-trees F a ^i where only 
one interaction a around i is retained, in formula one defines the effective magnetic 
field Hi acting on i, 

(o-i) = tanh( ) 8F i ) , (0.15) 
and deduces the latter from the effective fields in the sub-trees -Fo-w, 



Hi = 



E^ 



(0.16) 



This relation is represented schematically on the left panel of Fig. 0.5. These in turns 
can be computed from the effective fields in the sub-trees of the form Fj^. a , in which 
all interactions around the variable j except a is retained, 



— fa^i{{hj^ a }jeda\i) 



E 

b£di\a 



Ub- 



(0.17) 
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see the right panel of Fig. 0.5 for an illustration. The function f a -yi depends on the 
precise form of the interaction energy E a ; in the case of the p-spin interaction defined 
in Eq. (0.14) one finds 



1 

— arctanh 



tanh(/3J a ) tanh(^/i^ a ) 



(0.18) 



the present discussion applies to any kind of interactions E a , provided one takes the 
function / a _>i corresponding to the model under study. 

To summarize, the equilibrium properties of a tree factor graph model can be 
obtained by solving the recursion equations (0.17) for the effective fields {hi^ a , u a -n} 
living on the directed edges of the factor graph. These equations admit a unique 
solution which can be easily obtained by starting from the leaves of the tree (the 
variables with a unique neighboring interaction), where the boundary condition is 

hi-y a = 0. 

We can now proceed with the computations of the conditional averages (. . -)a , 
where the configurations of the variables in a subset B of the variables is fixed to 
g_ B . A moment of thought reveals that this is a simple modification of the previous 
computation. In fact one has to look for a set of effective fields {h~^ a , urL}, which 
depends on the imposed configuration g_ B . They obey the same recursive equations 
(0.17), but with an additional boundary condition: h~J^ a = Oi x oo if % G B. Indeed 
a fixed variable feels an infinite effective field (positive or negative according to the 
value it is fixed to). Once the fields {h~^ a , wff^} have been determined on all edges 
of the factor graph one obtains in a simple way the conditional magnetizations, for 
instance 

{a[) a _ B = tanh(/3iff B ) , Hf* = £ uf^ . (0.19) 

To complete the computation of the correlation function defined in equation (0.3), 
one has to average the correlation <Ji{a' i ) (jB with respect to Oi and a B , which are both 
parts of a reference equilibrium configuration a. The generation of an equilibrium 
configuration on a tree factor graph is most naturally done in a recursive, broad- 
casting like, way. This means that one has first to choose the value of a reference 
variable, say i, according to its marginal probability law. The latter is known on a 
tree, since the probability of Oi to be +1 in an equilibrium configuration is by defi- 
nition e@ Hi / (2 cosh(f3Hi)) = (1 + (<7j))/2. Once the value of <Tj is fixed in this way, 
the problem of the generation of the other variables decouples between the various 
branches emanating from i, thanks again to the tree structure of the factor graph. 
One can treat independently each interaction a around i, and generate the value of 
the other variables appearing in a, g_Q a \i, conditional on the value of en. The correct 
probability law for this generation reads 



M(£e \iki) = - ex P 



j—taGj 

j£da\i 



(0.20) 
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where z stands for a normalization constant. Once all the variables at distance 1 from 
i have been fixed, one can proceed in the same way with variables at distance 2, and 
so on and so forth until the complete configuration a has been generated. 

We can now wrap together the recursive computation of the conditional magnetiza- 
tions and the recursive broadcast generation of the reference equilibrium configuration. 
We shall indeed define Pi^aih'^a) as the (density of) probability that hfj^ a is equal 
to h' i _^ a when a B is part of an equilibrium configuration a drawn conditional on the 
value of (7j. A similar definition holds for QaU-i( u 'a^ri)> one can derive recursion rela- 
tions between these probability distributions, 

b£di\a \ b£di\a J 

QZUiK^i) = E ^GiaavkO / II dp Z«( h U«) (°- 22 ) 

Za a \i j£da\i 

complemented with the boundary condition Pj°^. a (/^_^ ) = S(h^_ >a — Cioo) when i G B. 
From the solution of these recursion equations one finally obtains 

°^ B) =^2^H~) J A «^U)tanh L £ 

-tanh 2 (^) . (0.23) 

The computation of the point-to-set correlation function we just presented is ex- 
act for any model whose factor graph is a finite tree. This is not precisely the case 
of problems on random graphs, as for instance the diluted p-spin model defined by 
Eqs. (0.12,0.13). In such models the factor graph is only locally tree-like, loops do 
exist but most of them have a length which diverge in the thermodynamic limit. In 
consequence the properties of random graph models can be inferred from such com- 
putations on finite trees, provided one takes into account the presence of the loops 
as self-consistent boundary conditions. This kind of reasoning has been extensively 
developed in the recent years under the name of the cavity method (32; 29), which is 
equivalent to the replica method (33) yet more practical for diluted systems. Following 
this methodology one can perform the average of the point-to-set correlation function 
(0.23) with respect to the quenched disorder of the model. Without entering the de- 
tails, this leads to functional recursion equations of the form P^ +i = F[Pi], where P is 
a probability distribution over real numbers (the h' in (0.21)), and the recursion is over 
the radius £ of the ball on which the correlation is computed. Numerical algorithms 
to deal with such equations between probability distributions have been extensively 
used to solve the equations of the cavity method. The idea, known in this context as 
population dynamics (32), is to approximate a probability distribution by a sample of 
elements, with an empirical distribution as close as possible to the true distribution. 

We show in the left panel of Fig. 0.6 the results of such a computation for the av- 
erage point-to-set correlation, as a function of the radius of the ball, in the case of the 
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p = 3 diluted model. The density parameter a = 1 is larger than the zero-temperature 
limit ad ~ 0.818 of the dynamic transition line (recall the sketch of Fig. 0.4). Indeed 
one clearly sees a plateau developing a diverging length upon decreasing the tem- 
perature towards Id (a) > 0. As discussed above the presence of this plateau makes 
the correlation length £(e) asymptotically independent of the precise choice of e, as 
long as it is smaller than the plateau value. This length is plotted in the right panel 
of Fig. 0.6 as a function of temperature. As the dashed line shows the divergence of 
the length occurs with a critical exponent 1/2, i.e. £(e) ~ (T — Td(a)) -1 / 2 . One can 
expect this exponent to be generic for almost all random graph models with such a 
discontinuous transition. It is indeed a natural consequence of the recursive structure 
of the equation underlying the computation, namely Pi+\ = F\Pg\. In the liquid phase 
the only fixed point of the equation P = F[P], reached when £ — ¥ oo, corresponds to 
a vanishing correlation. When the parameters (a, T) are taken to the glassy phase a 
bifurcation in this fixed point equation occurs, a non-trivial solution with a strictly 
positive correlation appears discontinuously. To simplify the discussion let us consider 
a toy model where the recursion is upon a real x instead of a function P, Xf+i = f(xe) 
(this corresponds actually to the zero temperature limit of the diluted p-spin model, 
see (34) for more details). The behaviour of the function / upon approaching such a 
discontinuous bifurcation is shown in Fig. 0.7. It is then a simple exercise to expand 
the function / when its parameter is slightly in the liquid phase, and its argument x is 
around the value of the close-by bifurcation, and to show in such a way that the length 
of the plateau in xi does indeed diverge with an exponent 1/2. The functional nature 
of the recursion in the general computation should not spoil this argument, one can 
thus expect to find gencrically such an exponent for the divergence of the point-to-set 
correlation length in mean-field diluted models. 

An important point should be emphasized here : for any kind of transition studied 
on a Bcthc lattice where a correlation length (shortest path on the tree) diverges as 
g-^tioc ( e o controlling the approach to the critical point), the mean-field exponent 
of the transition should be deduced as u m f = and not ^trec- Indeed the finite 
connectivity Bethe lattice should be pictured as embedded in an hypercubic lattice 
of very large dimension d. In this limit a path of length £ along the tree corresponds 
to a random walk (almost independently of the other branches of the tree) on the 
embedding lattice, with an end-to-end separation of order \f~£ in Euclidean distance 2 . 
The mean-field value of v for the transition discussed here is hence 1/4, we shall come 
back on this point in Sec. 0.4.1. 

0.3.4 Numerical simulations of diluted mean-field models 

The mean-field diluted models provide an ideal playground to assess the relevance of 
the bounds (0.8) between correlation lengths and correlation times: we have seen above 
how to compute the correlation length in the liquid phase and its divergence at the 
dynamical transition temperature T^. In addition the dynamics of these models can be 
easily investigated with Monte Carlo simulations (for out-of-cquilibrium dynamics see 
for instance (37; 38)). The equilibrium correlation times have been thus determined 

2 The validity of this argument, which is found for instance in (35; 36), can be checked on the 
simple example of a particle freely diffusing on a Bethe lattice. 
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Fig. 0.6 Left: Point-to-set correlation function for the p = 3 diluted mean-field model with 
a = 1, from left to right, T = 0.7, 0.62, 0.59, 0.58 and 0.575. Right: the point-to-set correlation 
length with e — 0.1, the dashed line is a fit to this data according to the asymptotic expression 
£ t (e) — C (T — T d (a))- 1/2 , with T d (a = 1) = 0.51695. 




Fig. 0.7 A representation of the bifurcation for the equation xi + i = f(xi). From bottom to 
top the parameters in / are in the liquid phase, precisely at the transition, in the glass phase. 

in (10), for a variation of the diluted p-spin model with constant degrees for the 
spins, and allowed to justify the claim of almost optimality of (0.8) without additional 
hypotheses. Indeed for T > Td one finds that the equilibrium correlation time remains 
finite in the thermodynamic limit, and diverges algebraically when T — > . Hence 
the lower bound of (0.8) could only been improved to include a dynamical exponent 
z and be promoted to r > £ z , with probably z > 2 for a large class of models. On the 
contrary in the low temperature phase, T < Td, it is the upper bound of (0.8) which 
is saturated (up to an improvement of the constant C2): the equilibrium correlation 
time diverges exponentially with the size N of the system, and the correlation volume 
is itself proportional to N. 

Coming back to the diluted p-spin model with Poissonian degrees, let us mention a 
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particular feature which is apparent in the sketch of Fig. 0.4. The dynamic transition 
line Td(a) reaches the zero temperature axis at a finite value ay- The neighborhood 
of the point (a = a<{,T = 0) thus exhibits a rich crossover phenomenon between 
low temperature activated dynamics controlled by energy barriers and positive tem- 
perature transitions of the schematic Mode Coupling Theory kind (which is an exact 
description of the a^co limit (39)). We refer the reader to (34) for more details. 

0.3.5 Connection with the replica symmetry breaking formalism 

We shall close this section on mean-field diluted systems with a discussion of the 
relationship between the computation of the point-to-set correlation function and the 
usual cavity formalism at the level of one step of replica symmetry breaking (1RSB), 
first unveiled in (40). The existence of such a connection should be expected by the 
reader accustomed to the fully-connected p-spin model: we argued that in the phase 
diagram sketched in Fig. 0.4 the transition line T<j(a), defined as the point where 
the point-to-set correlation length diverges, can be meaningfully called the dynamic 
transition. Indeed the lower bound in Eq. (0.8) implies a divergence of the equilibrium 
correlation time on this line. However the large a limit of the diluted p-spin model 
gives back the fully-connected one. For the latter the dynamic transition has been 
studied since a long time, and its appearance is known to be related to the static 
properties of the model studied by various methods, in particular the replica and the 
TAP ones. One of these characterizations is the existence of a non-trivial solution of 
the 1RSB equations with the Parisi breaking parameter m = 1. 

One can indeed recover this equivalent definition of the dynamic transition in 
generic mean-field diluted systems. Let us come back on the computation on a finite 
tree presented in Sec. 0.3.3. We introduced there the distribution P^\ a (h' i ^ a ) as the 
probability of hfj^ a with respect to the choice of the boundary configuration g_ B con- 
ditioned on <Ji . One can instead define the unconditional version of this distribution, 
Pi_ >a (/i^_ >a ), which is related to the conditional one by 



A similar definition holds for Qa-rtfa'a-ti)- A computation whose details can be found 
in (40) shows that the conditional equation (0.21) translates into 



This is nothing but the usual 1RSB equation with m = 1 (32) (in general the last term 
in parenthesis, called reweighting factor or free-energy shift, is raised to the power m). 
A similar observation holds for the equation on Q a ->i obtained from (0.23). The 1RSB 




(0.24) 
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equations are usually solved on a random factor graph without any explicit boundary 
B. In consequence the existence of a non-trivial solution of the usual 1RSB equations 
with m = 1 is equivalent to the presence of long-range point-to-set correlations (with 
a fictitious boundary sent to infinity). 

This computation provides also an alternative interpretation of the 1RSB equations 
for mean-field diluted systems. Pi_> a is usually considered as the distribution of an 
effective field with respect to the choice of a pure state of the system. This can be 
intuitively defined for a random graph of large but finite size as a portion of the 
configuration space sufficiently disconnected from the other pure states, yet a clear- 
cut mathematical definition is not easy to handle. One can alternatively interpret Pi^ a 
as a distribution over the boundary conditions of a tree of very large depth. This is 
precisely the computation done for the point-to-set correlation, i.e. the 1RSB equations 
with m = 1, where the boundary condition is drawn itself from the Gibbs-Boltzmann 
distribution. One can moreover generalize this interpretation to arbitrary values of m, 
biasing accordingly the probability distribution of the boundary conditions (41). 

0.4 Kac models 

A different class of systems in which the computation of PS correlation is possible 
is provided by Kac models. Kac models are classical tools of statistical physics to 
study the relationship between finite dimensional systems and their Mean Field coun- 
terparts (42; 43; 44). In these models one considers variables interacting through a 
potential of growing range ro and intensity decreasing in such a way as to keep fi- 
nite the total interaction strength of one particle with the surrounding environment. 
For any finite values of ro the general properties of finite dimensional system such 
as convexity of free-energy and absence of infinite life metastable states hold. There 
properties break down in mean-field thocry, which holds exactly in the limit of weak 
long range forces in which ro scales as the linear system size L and both lengths tend 
to infinity together. The relation between the two regimes can be understood studying 
the so-called Kac limit (42), where one considers a large interaction range ro — > oo 
but only after having taken the thermodynamic limit L — s- oo. Remarkably one finds 
that while the mean-field predictions hold in single phase regions, some of the typical 
pathologies of mean-field theory, such as the non convex free-energies are removed 
and the Maxwell construction emerges naturally. The resulting theory is an improved 
mean field theory that incorporarate the spatial dimension, with the possibility of 
describing inhomogeneous configurations and interfaces. For some systems (see, for 
instance, Refs. (45; 46)) it has been possible to devise asymptotic expansions around 
ro = oo which have led to the analysis of decay of metastable states or proofs of the 
existence of phase transitions. It is natural in the context of glassy systems to address 
the study of correlation lengths in the Kac limit. This provides a solvable limit case 
to confront with phenomcnological theories based on mean-field idcas(l; 9) and gives 
rise to a detailed picture with two relevant length scales. 

0.4.1 Definition 

Disordered models with Kac interactions were first described in (47), the case of spin 
glass models with pair interactions. Following the general philosophy of this chapter, 
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here we are interested in more general models with p-body interactions providing a 
minimal finite dimension modification of the fully connected p-spin model (0.12). 

We consider then spins er^ in the d dimensional hypercube i G A = {1, L} d with 
periodic boundary condition. The variables Oi can be of various nature, for example 
they can be Ising variables ui = ±1 or real variables subject to some local spherical 
constraint as we will specify below. The Kac p-spin Hamiltonian is defined by (48) 

E(a) = - J2 Ji 1 ,...,^i 1 ...<r ir (0.26) 

ii,i-2,...i p £k 

where, as for the mean-field model, the coupling constants Ji u ...,i are independent 
centered Gaussian variables, but now differently from the mean field model, the vari- 
ance depends on the mutual distances between the indexes according to: 

The function ip(x) can be chosen somehow arbitrarily, provided it fulfills the following 
constraints: 

• it is positive, il>(x) > 0. 

• it becomes negligible when its argument is much larger than 1. 

• it is normalized such that J d d x tp(x) = 1. 

These conditions guarantee that only spins within distances of order ro have an effec- 
tive direct interaction. The function tp can also be used to define a locally spherical 
model that without changing the physics from the Ising case, allows important simpli- 
fications in analytic studies. This is done considering Uj € R subject to the conditions 
^SjeA V' ( j a i = 1 m au points of space i. Note that the disorder distribu- 
tion implies that the Hamiltonian itself has a Gaussian distribution with covariance 
specified by 

E[E(a^)E(a^)]=J2Qi(s. W ^ (2) ) P > (0-28) 

i 

where Q%(g^ l \q^) — \ V\- ip ( ^— ^ ] aycr[ 2 ^ measures the similarity of the two 

configurations ct^ 1 -* and on a ball of scale ro around i. Models that compound 
different values of p can be defined considering gaussian Hamiltonian with correlations 
E[E(aW)E(aW)] = ^ <K<9i(2 (1) ,2 (2) )) ■ The function <j)(q) should satisfy suitable 
conditions to define a non-negative covariance matrix. 

Following the original suggestion by Kac (42), the scaling of the interaction with 
ro is chosen in such a way to insure the existence of the limit of the free-energy for 
ro — > oo. Different ways of taking this limit allows to study the relation between 
systems with finite range interactions and their mean field counterpart. If ro is chosen 
to scale with the system size ^ = L > 1 we have back the mean-field limit. All 
groups of p spins interact and the lattice structure and the space dimension become 
completely inessential. But one can consider a different procedure, called Kac limit, 
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where ro — > oo only after the thermodynamic limit L — > oo. A requirement for Mean 
Field to have a qualitative relevance in the description of finite range systems is a 
smooth ro cross-over from ro <C L to rp ~ L (L 3> 1). Rigorous results connecting 
the Kac limit to the mean- field limit have been proved in (50; 48), under the technical 
hypothesis of parity of the function <j), cf>(q) = <p(—q). These findings can be informally 
and qualitatively stated as follows: 

• The free-energy of the Kac model tends for all temperature to the free-energy of 
the corresponding mean-field model in the Kac limit. 

• In the same limit, the structure of the correlation functions on length scales of 
order of ro approaches the structure of mean-field correlations. 

The first point is a continuity property, it tells that the mean field limit is a good 
starting point in the evaluation of the free-energy for large ro. Of course this does 
nor necessarily imply that the phase transitions found in mean field are also present 
for some finite ro. For example, it is clear that in ID the free-energy is an analytic 
function of temperature for all finite ro, despite the fact that the large ro limiting 
function exhibit a phase transition. 

The second point clarifies that mean-field approximation is at least capable to 
capture the local properties of the system. This is particularly interesting at low tem- 
perature where it implies that the system behaves as an ideal glass at least locally 
on length scales of order ro- The possibility that this description extends to scales 
much larger than ro in some finite space dimension remains an open question. To go 
beyond these general rigorous results and study quantitatively the behavior of corre- 
lation lengths, tools of theoretical physics and non-rigorous techniques are needed as 
it will be explained in the next sections. 

0.4.2 Effective potential: a Landau free-energy functional for glasses 

In this section we would like to review the basic tool of the glassy effective poten- 
tial (52), which has been instrumental in the computation of point-to-sct correlation 
lengths in the Kac model (58). 

The effective potential is a Landau kind of free-energy as function of an order pa- 
rameter explicitly devised to be sensitive to stable or metastable glassy states. While 
this construction is in principle very general, for sake of simplicity we first expose it in 
the case of fully connected models, and later explain how to extend the analysis to Kac 
models. The basic idea is that metastable states can be studied through restrictions of 
the Boltzmann-Gibbs measure to the surroundings of a representative configuration. 
One then observes that in glassy systems where equilibrium is composed by an ex- 
tensive number of metastable states, any equilibrium configuration can be taken as a 
representative configuration. Let us then choose a random configuration a according 
to the equilibrium measure, and consider a second system a which is in a constrained 
equilibrium where the overlap with the reference configuration q(g_(°\q_) = ^ fjof 
is fixed to some preassigned values p. 3 

3 Note that we used already the letter p to define the p-spin Hamiltonian. It will be clear from the 
context to which quantity the notation refers. 
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The free-energy of this second system, for fixed a<°> is 

F[P, ^ (0) ] = -§ Z{p,<L {0) ] (0.29) 
2b,2 (0) ]=E e "^II%fe (0) ^)-f)- (0-30) 

In the thermodynamic limit this is is self-averaging with respect to the quenched 
disorder in the interactions and the choice of the equilibrium configuration (7°, 

F\p]=E(F\p,c L M}) , (0.31) 

where by (•) we denoted the average with respect to the equilibrium distribution 
for cr^ -* and by E the average with respect to the quenched disorder. The function 
F[p] can be used as a Landau free-energy as a function of the order parameter to 
detect and study glassy states. It is customary to define the effective potential V\p] 
subtracting the unconstrained value of the free-energy, V[p] = F[p] — F. This is the 
large deviation function allowing to evaluate the probability F[p] of the overlap with 
a reference configuration 

P[p] = exp(-pv\p]) . (0.32) 

The analysis of the effective potential, which requires the computation of the double 
average in (0.31) can be performed in mean-field theory through the replica method, 
or alternatively the cavity method. For the family of models we are interested in, the 
most general ansatz needed to describe V[p] is of a 1RSB form, that accounts for the 
possibility of having multiple metastable states. In fully connected models this takes 
a simpler form with respect to the diluted models described in section 0.3.5. This is 
just parametrized by the Parisi breaking parameter m and two real variables q\ and 
go, which represent respectively the typical value of the overlap between configuration 
in the same metastable state and in different ones (2). The parameter m can be tuned 
to select different families of metastable states with different free-energies. If the equi- 
librium phase is composed by an extensive number of metastable state, corresponding 
to a non-zero configurational entropy m takes the value m = 1. 

In the spherical p-spin model, defined by the Hamiltonian (0.12), the effective 
potential can be computed in a closed form as a saddle point over qi, q^: 

V\p]=max.V(p,qi,qo,m) (0.33) 

9l,9o 

V(p,qi,qo,m) = --/3[2^(p) - mc/)(q ) - (1 - m)4>{qi)] + 

^ f i — 7T~ — r 2 + ~ — ~ lo ^ 1 ~ 21] _ ~ ^st 1 - 0- - m )<?i - m 1o] 

Ip \ 1 — (1 — m)qx — mqo m m 

where we denoted 4>(q) = ^q p . The shape of the function V\p], which depends on 
temperature, is sensitive to the existence of metastable states. Indeed, as can be seen 
in Fig. 0.8, one can identify three model dependent temperatures of interest where 
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qualitative changes take place. Two of them are the mean-field dynamical transition 
temperature Td and the static transition temperature T c . In addition there is a third 
temperature T*, with T* > T d > T c , first identified in (52). Above T* the function 
V is a convex function with a single minimum at p = 0. At T* and inflection point 
appears, and below that temperature the potential is non-convex. For temperature 
between Td and T* , the function continues to have a single minimum for p = 0. At T d 
a local minimum at a value p = qEA > develops. In the interval [T c ,Ti], the point 
p = is still the absolute minimum of V. The two minima structure below Td reflects 
the partition of the equilibrium measure in disjoint metastable states. The value qEA is 
the typical overlap between configurations belonging to the same metastable state. For 
P = qEA 3_ is in the state specified by . Different metastable states have zero mutual 
overlap. For p = all but the metastable state specified by contribute to the free- 
energy and V(0) = 0. Correspondingly, the difference in free-energy between the two 
minima equals the system's configurational entropy Eoo(T) multiplied by temperature. 
The configurational entropy vanishes linearly on approaching T c , Eoo(T) ~ T — T c and 
the two minima become degenerate. Below that temperature the mean field model is 
in an ideal glassy state and the two minima remain degenerate. 




Overlap p 

Fig. 0.8 The function V(p) at different temperatures. For comparison with the case 
of Kac model we consider a Hamiltonian with two body and 4 body interactions with 
<j>(p) = 1/2(0.1 x p 2 +p 4 ). From top to bottom T = 0.703486 > T* , T = T* = 0.633137, 
T — Td — 0.57525, T = 0.558049, T = T c = 0.541847. The function is convex for T > T* . 
It has an inflection point with positive slope for Td < T < T*. In the interval T < T < Td, 
V(q) has a local minimum for a temperature dependent value p = qEA- The difference 
V{q E A) - V(0) is (T times) the bulk configurational entropy V{q E A) - V(0) = TE oc (T). 

We now generalize the previous construction to Kac models. As before, we use 
the overlap with an equilibrium reference configuration as an order parameter. 
In this case however, we are interested in considering the free-energy as a functional 
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of an order parameter profile that is space dependent. In view of considering large 
values of the interaction range rg it is natural to coarse-grain the order parameter on 
scales 5 such that S ~ dx r<j ^> 1 with dx <C 1 (59; 60). This allows to rescale all 
lengths x — i/r D which are then measured in units of the interaction range rg. Let us 
denote q x (s.^°\3.) = jx J2i^B <T i°'' cr i the local overlap in a square cell B x of linear size 
5 labeled by x between two configurations cr*- -* and a. In complete analogy with the 
long-range case, one can define an effective potential functional W^[p(x)] = F\p(x)] — F 
as the average free-energy excess for imposing an overlap profile p(x). 

FW)\ = -pElog Z\p(x),aW] (0.34) 
Z[p{x) : a^] = Y j e- pE ^X[5{ qx (^\c L )-p{x)). (0.35) 

q_ x 

We can now relate the point-to-set correlation function to the effective potential. Con- 
sider the correlation function between point and the set B considering profiles such 
that p(x) = 1 is fixed for x £ B. Observing that the probability of a profile can be 
written as P[p(x)] = exp(— W(p(x))) wc find 

«P(0)> («>) = {P r (x) ^ 1;XgJ?} S \ RW < „ ■ (0-36) 
/ w=liJ6B} ^)e-W] 

For general values of ro the evaluation of the large deviation functional W[p(x)] is 
a formidable task. Fortunately, in the Kac limit ro — > oo the computation becomes 
simpler, it can be stated as a variational problem that generalize to spatially varying 
functions the saddle point evaluation of (0.34). The region of temperatures T > T c can 
be analyzed through the 1RSB theory that we have discussed in several contexts in this 
chapter. This generalizes the description of the fully-connected case in the fact that the 
overlap variational parameters become position dependent functions qi(x) and qo(x) 
(the parameter m remains space independent). We give the explicit expression of the 
effective potential for the p-spin model (0.26) with local spherical constraint, where 
the dependence on the various parameters can be written in a closed form 

W[p(x)\ = max r d [ d d x (-U{2[4>(p*{x)) ~ *(p(*))] " "#(<£(*)) - 0(«b(»))] 
-(l-mMqt(x))-ct>(q 1 (x))}}+V(p(x),q 1 (x),qo(x),m)) , (0.37) 

We have introduced the notation f*(x) = J d d y ip(\x — y \)f(y) for the convolution 
with t/; and the function V(p,qi,qo,rn) is given by the mean field expression (0.34). 
Mutating the terminology from mechanics we will call action the argument of the 
max in (0.37). We notice that for profiles that are inhomogeneous on limited spatial 
extension W is proportional to the interaction volume Tq. The expression (0.37) can 
be simplified using a (physically harmless) lowest order gradient expansion for the 
convolution <j>(f*(x)) = <j>{f(x)) + eft {j{x))V 2 f{x) with c = ±Jd d z ip(\z\)z 2 . This 
approximation does not affect the physics of the problem while the action takes the 
familiar field-theoretical form of the integral of a Lagrangian density, dependent on 
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the various overlap parameters and their spatial derivatives. Notice however that for 
the expansion to be well behaved one needs 4>"(q) > for all q. This is not verified for 
pure p-spin interactions if p > 3 which has 4>"(0) = 0. For this reason the numerical 
computations based on this truncation have used functions <f> of the kind <p(q) = 
\{q A + aq 2 ), which correspond to a mixed p-spin model with two body and four body 
interactions. For small enough values of a the model has the same phenomenology as 
the pure p-spin. Figures are presented for the value a = 0.1. 

Thanks to the proportionality of the action to the interaction volume the func- 
tional integral in the expression (0.36) of the PS correlations can be evaluated by 
saddle points. We will be then interested in profiles p(x) that are stationary points of 
the action = for x £ B with boundary condition p(x) = 1 for x £ B. These in- 
clude the absolute minima, that represent equilibrium profiles and directly determine 
the PS correlation function, but also relative minima, saddle points and maxima that 
we will associate to metastable states and barriers. 

The condition of equilibrium (0.6) discussed in section 0.2.2 has interesting impli- 
cations on the mutual relation between the overlap functions q\{x), qo(x) and p{x). 
Consider the local overlap q x (£.',^'), between two replicas with the same constraint 
er' ) on B, and its fc-th moment E(q x (oi',g_") k )- it is possible to show (2) that within 
the 1RSB formalism 

lim E(q x (a' lS !') k ) = (1 - m)q k (x) + mq%(x) (0.38) 

r -»oo 

where qi(x) and qo(x) are the solution of the maximization problem (0.37). On the 
other hand, according to (0.6), this expectation value should be simply equal to 
\im. r0 ^ oo K(q :r (a_ < - a \a_') k ) — p{x) k , thus leading for all k and x to the identity (1 — 
m)q k (x) + m<7g (x) = p{x) k . Two different ways of satisfying this identities are of phys- 
ical relevance: 1) a simple solution describing equilibrium in a single metastable state, 
such that for all x 

qi (x) = q (x) = p(x) . (0.39) 
2) A solution with m — > 1, for which 

3o (x) =p(x) , 

qi (x)>p(x). (0.40) 

This solution (see below) is relevant when ergodicity is broken and the constrained 
equilibrium is composed by multiple metastable states with a non-zero configurational 
entropy S. According to general principles (62) the configurational entropy can be 
computed from a solution with the form (0.40) as S = — P^^^\ m =i- 

Non-equilibrium stationary profiles p(x), associated to transition states and bar- 
riers, will be relevant in our analysis. For these non-equilibrium solutions, there are 
no restrictions of principle on the functions qi(x) and qo(x). Of course the condition 
qo(x) = qi(x) = p(x) can still be used as an approximation. Despite the fact that in 
some regions this gives rise to clear physical inconsistencies we believe that it predicts 
the correct scaling of the free-energy of metastable and barrier states with temperature 
and interaction range. More complex solutions, studied in (69) confirm this point. 
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Let us stick to the forms (0.39) and (0.40). In both mentioned cases, one gets 
great simplifications in the variational problem one needs to solve. For example for 
the ergodic equilibrium profiles (0.39) the effective potential, using the mentioned 
gradient expansion for the convolution reads: 



In the 1RSB solution (0.40) an additional qi dependent term proportional to 1 — m 
appears that allow to compute the configurational entropy. It is natural to look for 
solutions respecting the symmetries imposed by the set B. We will discuss mainly 
the spherical geometry, where the set B is just the exterior of a spherical cavity, 
B = {x : \x\ > £} and solutions should only depend on \x\. Another interesting 
choice is the one of planar geometry, where B = {x : \xi\ > £}, with x\ the first 
component of the d dimensional vector x, solutions should only depend on x±. In both 
cases it is interesting to study the behavior of the solutions as a function of the cavity 
size £, that, remember, we measure in units of ro as all the lengths. Physical saddle 
points should be such that in both cases the profile is flat in the center of the cavity 
Vp(0) = 0. In both geometries the minimum equations for W can be easily numerically 
integrated, yielding solutions qualitatively independent from the physical dimension 
in the spherical case, and strictly independent of the physical dimension in the planar 
case. The main qualitative features can be obtained in the planar case where the 
overlap profile can be obtained by quadrature. In fact this case is formally identical 
to the problem of a ID Newtonian particle with coordinate p evolving in time in a 
potential equal to — V(p). Unless otherwise stated the results we will discuss, and in 
particular the behavior of the different lengths that will be identified, are independent 
of the geometry and the dimension. We will denote p(x; £) the solution p(x) of the 
variational problem for the boundary conditions fixed at distance I. 

0.4.3 The physical picture 

The nature of the solutions of the stationary point of (0.41) and the behavior of the 
PS correlations as a function of the cavity size £ reflect the shape of the mean-field 
potential V{p) that we have discussed in the previous section. 

In the "very high" temperature region T > T* for each value of £, both in the 
spherical and planar geometries, there is a unique solution of the field equations with 
the imposed boundary conditions. As a function of £ the PS function p(0; £) decreases 
continually from 1 to 0. 

At T = T* an instability develops. Below that temperature one finds a temperature 
dependent interval of lengths / = [^o,^d] where the field equations admit three distinct 
solutions. For temperatures in the interval (Td, T*) both lengths are finite, while £& = 
oo for T < Td. The behavior of p(0;£) for different temperatures is depicted in the 
figure 0.9. The figure also shows the complete profile p[x\£) as a function of x for 
specific values of £ and the temperature in the three solutions. According to the value 
of p(0;£), we call the three profiles Low Overlap Solution (LOS), Medium Overlap 




V(p) = - - W) + log(l - P(*))] • 



(0.41) 
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Solution (MOS) and High Overlap Solution (HOS). These have different free-energies. 
The HOS and the LOS one correspond to free-energy absolute or relative minima, the 
MOS one is a maximum with respect to p(0; £). 

The structure of the solutions allows to identify three relevant lengths which grow 
for decreasing temperatures: two metastability lengths £q and £d defining the interval 
in which three solutions exist, and a thermodynamic length £ c , with £q < £ c < £&, for 
which the HOS and the LOS have the same free-energy. The presence of two free-energy 
minima separated by a barrier naturally interpreted a metastability phenomenon. The 
typical relaxation time depends on size of the cavity i and diverges exponentially for 
r -> oo if £ < t < id- 

We divide the remaining of this section in the discussion of the equilibrium solu- 
tions, the metastable solutions and the unstable solutions for T < T*. 

Equilibrium and metastable solutions. Let us consider the behavior of the system in 
spherical geometry as a function of the ball radius I. Below the length £q the HOS is 
the only minimizcr of the free-energy functional, there are no metastable states poorly 
correlated with the boundary conditions. This absence has consequences in dynamics: 
any configuration in B chosen as initial condition for the dynamics, evolves in an 
tq independent time to the HOS. Above this length, in the interval £ € {io,£ c ) the 
LOS makes its appearance as a free-energy local minimum. Though the HOS is still 
thermodynamically dominant, metastable states exist such that if they are chosen as 
initial conditions in relaxation dynamics, they will take a time exponentially large in 
the interaction volume r$ to relax to equilibrium. 

The difference in free-energy between the HOS and the LOS decreases upon increas- 
ing £, until, for £ — £ c the two free-energy minima are degenerate. As we will discuss 
later, the length £ c is growing for decreasing temperatures, and diverges for T — > T c . 
Above this length the LOS becomes the dominating thermodynamic state. There is a 
phase transition in £ separating a phase strongly correlated with the boundary below 
£ c , from an essentially uncorrelated phase above. The PS correlation function which 
should be identified with the value of p(0; £) in the thermodynamically dominant solu- 
tion exhibits a first order jump in £ = £ c . Below £ c the system is thermodynamically in 
an ideal glassy amorphous state where, despite no sign of intrinsic order, the boundary 
conditions determine the bulk configuration. Conversely, above £ c , the system is in a 
liquid situation where the boundary has little influence on the bulk. The presence of 
the transition is reminiscent of the analysis of PS function by Biroli and Bouchaud (9), 
with £ c thus identified with the mosaic length of phenomcnological RFOT (1). This 
theory postulates a competition between interfacial energy and configurational en- 
tropy, the contributions compensating exactly at the transition. However, as we have 
discussed for generic observables in sec. 0.2.2 (see formula (0.4)), the average energy 
density for the constrained system is equal in all points of space to its unconstrained 
equilibrium value. This implies in particular that at the coexistence point both phases 
should have the same energy density, and therefore entropy, in each point of space. The 
transition should therefore depend on a different mechanism. Further analysis of the 
HOS and the LOS presented below reveals that it has the character of a configurational 
entropy crisis similar to the one that happens at T c for the infinite system. 
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In the region £ £ (i c ,^d), the HOS appears to be a local minimum of the free- 
energy functional. The length £d marks the limit of metastability of the HOS and was 
called dynamical length in (58). For r$ — ► oo and £ £ (l c ,£d)i a system prepared in 
the configuration at an initial time t = would not be able to dynamically relax 
towards the thermodynamic equilibrium profile. For large but finite r$ the relaxation 
time (for the dynamics in a ball of radius i) to the LOS profile can be expected to 
behave exponentially, t{tq,1, T) ~ exp(rp^(^, T)). In other words, for ro — > oo, while 
the relaxation time is divergent, the barrier (counted in units of the interaction volume) 
lim r . _ ! . OC) 4j logr(ro, £, T) = B(£,T) remains strictly positive. Conversely, for £ larger 
than £d the relaxation time would remain finite for ro — > oo. id appears as the minimal 
length for which relaxation is possible without barrier jumping. As Fig. 0.10 shows, 
the dynamic length increases monotonically and diverges as T decreases from T* on 
approaching Td. Generically in the class of models we are interested in, it behaves as 



at least below dimension 6. Interestingly, this scaling, first found in (1), coincides with 
the one obtained for the length £4 quantifying the typical extension of collective motion 
in MCT (7; 8) relaxation processes. This seems natural since, in the r — > 00 limit a 
set of ideal inhomogeneous MCT equations can be shown to describe the dynamics of 
the Kac model (63). A simple scaling argument accounts for the above behavior in the 
present context. For temperatures close to T d and dimension smaller than d c = 6, the 
properties of the solutions only depend on V in the vicinity of the inflection point (64). 
The potential can then be expanded to the cubic order around j>d, its inflection point 
at Td, &sV(p,T) = V(pd 7 T) + ea(p — pd) + b(p — pd) 3 , where generically the coefficient 
of the linear term is e = T — Td while a and b are positive model-dependent parameters 
with a weaker temperature dependence. Writing p — pd — <p this gives rise to a cubic 
field theory 



whose critical properties dominate the behavior of observable quantities for e — J- 0. In 
particular, equating the order of magnitude of the different terms in the Lagrangian 
etp ~ ip 3 ~ (V<y9) 2 immediately leads to a length £d ~ e -1 / 4 . It has been noticed in the 
context of spinodal theories that the expansion (0.43) cannot be used above dimension 
6, and the critical point has a non- universal character (64). 

The nature of the LOS is clarified considering the 1RSB solution for m = 1 (0.40) 
we mentioned in the previous section. In the region I > £d the LOS is the unique free- 
energy minimizer and it represents an ergodic state. Conversely, for £ c < £ < £d this so- 
lution describes a non-ergodic state composed by a collection of J\f ~ exp(ro^ d S(£, T)) 
metastable states. The free-energy difference between the HOS and the LOS, which, 
as previously observed, should be purely entropic, in fact coincides with (T times) the 
configurational entropy in the constrained system, denoted £ d T,(£,T). In other words 
there is an exponentially large number Af of metastable states compatible with the 
boundary condition. The HOS just represents one of them, on the same foot as all the 
others. The behavior of S(£, T) is depicted in the lower panel of Fig. 0.11, which shows 



*d~(!T-T d )-a , 



(0.42) 




(0.43) 
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curves growing with £, crossing the value for £ = £ c . If T > the curves stop at 
the dynamical length £^ beyond which the system recovers ergodicity. For T < the 
curves tend to the bulk value of the mean- field configurational entropy E^fT). The 
dependence of T) upon £ can in both regions be described in a good approximation 
as E(T,£) = A(T)(1 - £ c /£), where for T < T d , A(T) = Soo(T). Upon decreasing the 
values of £ towards £ c , the transition to the confined state is driven by the reduction 
to zero of the configurational entropy. This is decreased by a term that scales as £ 1 
relative to the total configurational entropy £ d T,(£, T). 

Let us mention without entering in a detailed discussion that in addition to the 
solutions we discussed, low-overlap RSB solutions with m ^ 1 (not respecting the 
conditional equilibrium conditions) were also found in (58). These represent families 
of mctastable states in the system with free-energy higher then the equilibrium one. 
For £ < £ c the solution with m = 1 gives a negative configurational entropy and 
is therefore inconsistent. The dominating metastable low overlap state is one of the 
RSB m/1 solutions just mentioned. It has zero configurational entropy and positive 
energy difference with the HOS. 




Coordinate Ixl 

Fig. 0.9 Upper panel: values of the PS function p(0, £) for three values of the temper- 
ature in the different solutions: T = 0.676 > T* = 0.633, T = 0.610, in the interval 
[T d ,T*] = [0.575,0.633] and T = 0.569, in the interval [T c ,T d ] = [0.541,0.575]. In the first 
case there is a unique solution to for all values of I. In the second case three solutions are 
present for I in the interval [£o,^d] = [17.31,31.57]. In the third case three solutions are 
present for I > £o = 19.47. The arrow mark the points £ — £ c . Below £ c the HOS is the 
thermodynamically favoured state and the LOS is metastable, above that length the roles 
are interchanged. For the T = 0.610, £ c = 20.4, for T = 0.569, 4 = 26.7. Lower panel: The 
overlap profiles p(x;£) corresponding to the the three solutions for £ — 20 and T = 0.610. 

The scenario for the various glass transitions as a function of £ in the LOS resem- 
bles the one of the Mean-Field or diluted p-spin models as a function of temperature. 
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Fig. 0.10 Thermodynamic (filled squares) and dynamic (empty squares) lengths £ c and 
id- The vertical lines correspond to the dynamic and thermodynamic glass transitions 
T c = 0.541847 and T d = 0.57525. The lengths behave respectively as 4 ~ (T - T)- 1 close 
to T c and as id ~ (T — Ta) -1 ' 4 close to Td. In the inset, the lengths as functions of (T — Td) 
(for id) and (T — T c ) (for i c ). The dotted lines have slope (respectively) 1/4 and 1. 

We have an ergodic phase for large £ > id, an intermediate phase with a finite con- 
figurational entropy for £ c < I < £d, and a zero configurational entropy 1RSB phase 
at low £ < £ c . We notice that the point at £d has the character of a typical dynamical 
transition point of 1RSB systems. This has specifc consequences on the dynamical 
properties (for the constrained dynamics in a ball of radius £) when £ gets close to £d, 
which should be similar to the one of the mean-field model close to Td. In particular, 
for ro — > oo, the relaxation time grows and diverges if the system size goes to £d from 
above. On the same foot one can draw consequences on the relaxation time within 
the metastable states dominating below £<\. One can argue that the stability of these 
states decreases upon increasing £ until they become marginal exactly at id (65). Ac- 
cordingly, the relaxation time within a metastable state diverges if £ tends to id from 
below. 

As already mentioned the dynamics of the model in the Kac limit is exactly de- 
scribed by inhomogencous MCT equations (63), so that the bulk relaxation time 
behaves as r(T) ~ e^ 1 close to Td, with a model dependent exponent 7 (66). For 
£ > id ~ e" 1 / 4 one can expect the scaling 

r(£,T) = ±-r(i/id), (0.44) 

where the scaling function f(u) tends to a constant at large arguments, and has a 
power law behavior f (u) ~ (u — 1)~ 5 for u > 1. The exponent 5 should be in principle 
derived from the dynamical MCT-like equations (63) and as 7 should depend on the 
choice of the model. An analogous scaling should be expected below £d with t(£, T) 
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interpreted as the relaxation time within a state, for the dynamics inside a ball of 
radius I. 



0.01 
0.001 
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Ball radius 1 



£ 0.006 
I 0.004 
J 0.002 

I o 

a -°- 002 

16 18 20 22 24 26 28 30 
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Fig. 0.11 Upper panel: The barrier density as a function of the ball radius £ for different 
temperatures. From bottom to top T = 0.5819, 0.5816, 0.5812, 0.5809, 0.5805, all larger then 
Td = 0.57525. Lower temperatures correspond to higher barriers. The barrier goes to zero 
for £ — > £d- The dotted line a allows the comparison with the power law £~ 6 expected from 
scaling close to Td. Lower panel: Configurational entropy density as a function of I for the same 
temperatures as in the upper panel. The higher curves correspond to the higher temperatures. 
The configurational entropy curves touch zero at the mosaic length £ c and terminate at the 
dynamical length £&. Both lengths increase for decreasing temperatures. Below the length £& 
the correct solution should include RSB. 

The thermodynamic mosaic length £ c , whose temperature behavior is shown in 
Fig. 0.10, displays no singularity at Td. This is coherent with the fact that for any 
finite ro the dynamical transition is roundcd-off and the relaxation time remains fi- 
nite, and is in agreement with the inequalities (0.8) of Sec. 0.2.3. Its behavior becomes 
singular near the ideal glass transition point T > T c . Here the two minima of V are 
almost degenerate, and 4 is large. The solutions of the field equations in spherical ge- 
ometry for d > 1 can be obtained through the "thin wall approximation" of nuclcation 
theory (67). This self-consistently neglects the thickness of the interface region, where 
the solution passes form values of p close to qEA to values close to zero, this thickness 
being small with respect to £ c . Let us consider £ > £ c . The transition is driven by a 
bulk contribution to the configurational entropy E oc (T)€ d and an almost temperature 
independent surface entropy reduction — Y l d ~ x . The value of Y can be computed in 
analogy with standard nuclcation theory as the action of the "instantonic solution" to 
the field equation that connects the two degenerate minima of V(p) at T = T c : 
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Y= Sd r dP - 2v ^ 



(0.45) 



T=T C 



where Sd is the surface of the unitary sphere in dimension d, Sd = 2it d / 2 T(d/2). The 
length £ c is then found from the cancellation of the total configuration entropy, 

£ d Z(£, T) = £oo(T)£ d - Yi d ~ l , (0.46) 

namely £ c = YfS^T). The length £ c diverges at T c , with the scaling 

Y Y 

(0.47) 



Soo(T) T-T c 

We will see in the next sub-section that this is also the typical spatial extension of 
barrier states close to T c . Notice that the negative contribution to (0.46 ) is propor- 
tional to the ball surface. As we discuss in sec. 0.4.3 phenomenological theories (1; 9) 
suggest that interface reduction terms could scale as £ e with 6 < d — 1. As all qual- 
itative results we present the value 9 = d — 1 can be checked to be independent on 
the choice of the model in the 1RSB class or the various approximations like the thin 
wall approximation or the gradient expansion of the interaction that we have used 
to simplify the analytic treatment. Our derivation, however, depends crucially on the 
Kac limit, that in our analysis preceeds the limit T — > T c . It is certainly possible that 
for finite ro a non trivial cross-over could change this exponent. 

The barrier state. Let us complete our analysis with a discussion of the medium 
overlap solution (MOS). In figures 0.9 the value of p(0; £) plotted in the MOS has been 
obtained supposing conditional equilibrium as an approximation. Unfortunately, it is 
difficult to go beyond this solution using numerical integration. Since by definition the 
MOS is an unstable saddle of W, iterative methods that work well for minima, simply 
do not find the solution. Experience with the replica method shows however that such 
replica symmetric solutions often give good approximations for the free-energy and 
other quantities even in regions where the exact description should include RSB. 

We analyze below the regimes of temperatures close to Td and T c where scaling 
can be expected and analytic arguments allow to go beyond numerical integration of 
the field equations. We believe that our approximate solution while possibly failing in 
numeric prefactors, gives back the correct scaling behavior describing the dependence 
of the MOS on system size and temperature for large ro. 

The MOS is interpreted physically as a dynamical barrier state to be overcome 
during relaxation (68). The equilibrium activation barrier B(£,T) can be estimated as 
the free-energy difference between the MOS and the HOS 

r d B{£, T) = Wmos - Whos ■ (0.48) 

This is, for all temperatures, a decreasing function of the system size £. For temper- 
atures in the range [Td,T*] the barrier vanishes at £<±, coherently activation is not 
needed for relaxation in large systems. For temperatures below Td activation is re- 
quired at all scales. The barrier for bulk relaxation can be obtained as the large £ limit 
of B(£,T). 
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Close to Td the behavior of the MOS and the HOS become independent of the 
details of the model. Since we do not expect RSB to affect the scaling, we can obtain 
the behavior of the various quantities through dimensional analysis of the cubic field 
theory of Eq. (0.43). The main results of this analysis are that 

• For T > Td the barrier is non zero for t < I d and admits the scaling form 

B(£,T) =£ d - 6 b+(£/£ d ) , (0.49) 

where the scaling function is such that the barrier vanishes linearly as £ — > £ c \: 
b+(x) sa Const, x (1 — x) for x — > 1 and b+(x) = for x > 1. The behavior of the 
barrier in this region is shown in the upper panel of fig. 0.11. 

• For T < Td the barrier admits the scaling form 

B{£,T) = \ef^b.{£\e\ 1 ' i ) , (0.50) 

where the function 6_ (x) is a decreasing function of x that admit a positive limit 
for x — > oo. As already noticed, from (0.50) we can obtain the barrier for relaxation 
in the bulk as the limit of large £ of B(£, T). This coincides with a direct estimate 
of the barrier in (61; 69) where boundary conditions are imposed at infinity and 
scales as Boo{T) ~ | e | 3 . The scaling (0.50) shows that the asymptotic value is 
reached on lengths of order £b ~ lei" 1 / 4 . This is then the typical spacial extension 
of dominating barrier modes governing the relaxation in the bulk, which becomes 
large on approaching Td from below. 

Let us now turn to temperatures close to the static transition temperature T c . 
Analogously to the case of the equilibrium solutions of the previous subsection, for 
small T — T c nucleation arguments can be used to study the barrier. For large £, the 
barrier solution p(x; £) has the shape of a bubble of small overlap values in a sea where 
q(x;£) w qEA- This is surrounded by a crown around the boundary of B where p 
decays from 1 to ps qEA- For d > 1 the typical bubble radius £b, and the barrier value 
can be just estimated from the point of maximum of expression (0.46) 4 . This gives, as 
announced, a common scaling for £^ and £ c , 

'b-^-Cr-r^, (0-51) 

and for the barrier the behavior 

B{£, T) ~ F d (S 00 (T))-^- 1 ' - (T - TJ-V-V , (0.52) 

which is divergent for d > 1 when T — > T c . Eq. (0.52) has the form of a modified 
Adams- Gibbs inverse proportionality relation between the relaxation barrier and the 
configurational entropy(70). In turn, using Soo(T) ~ T — T c , one finds a modified 
Vogel-Fulcher law 5 with an exponent that depends on dimensionality, with associated 
lower critical dimension dhCD = 1- We remark the coincidence of this result in the 
Kac limit with the naive analysis of (1). The comment we made at the end of section 
(0.4.3) on the exponent can be verbatim applied here to the modified Vogel-Fulcher 
exponent. 

4 The Id case requires a special analysis and one finds that £ h ~ — log(T — T c ). 
5 The usual Vogel-Fulcher form posits a barrier proportional to (T — Tc) -1 . 
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Correlations in wall geometry. Before concluding our discussion on the Kac limit, 
let's have a closer look to the three solutions for large £ ^> £ c and T >T C . The large 
£ limit is equivalent to imposing the boundary condition on a single planar wall. This 
condition has been object of recent numerical experiments (18) and deserves a short 
comment. 

The question here is up to which distance the boundary has an effect on the 
configuration of the system, i.e. what is the typical distance from the boundary over 
which p(x;£) decays to zero and if this distance coincides with £ c . The behavior of 
the three solutions in this regime can be easily worked out from the field equation. 
The metastable HOS never looses correlations: after a rapid decay from 1 on the wall, 
which is common with the MOS and the LOS, stabilizes to p(x; £) ~ qEA- The unstable 
MOS persists to values p(x; I) « qEA up to distances of order I from the wall. It decays 
to zero at a distance £ c from the center and is similar to a ball of radius £ c of small 
overlap in a sea of large overlap values. More interesting is the behavior of the LOS 
which gives the equilibrium profile. In this case the overlap decays to low values over 
distances of order of a "wetting" length £ w ~ — log(T — T c ) from the boundary before 
decaying to values close to zero. We conclude that the semi-infinite wall geometry fails 
to identify the correlations £ c while it identifies a much smaller length £ w . 

Beyond the Kac limit. The results presented so far are strictly valid in the Kac limit. 
We briefly discuss here some possible scenarios for ro large but finite. 

We would like first to discuss the rounding off of the dynamical (MCT) transition. 
The analysis can be performed in analogy with ferromagnetic models with Kac interac- 
tions below the lower critical dimension where "finite range scaling" (FRS) holds (71). 
The basic hypothesis of FRS is that the interaction range acts as a cut-off to critical 
behavior in a similar way as finite size does for usual critical phenomena. This idea has 
been employed in a Id Kac spin glass with a continuous transition in the Kac limit, to 
predict the rate of growth of the correlation length with ro at the mean-field critical 
temperature of the model (72). The ro — > oo dynamical critical point has a mean- field 
character, so that one can expect mean-field scaling. An observable O which depends 
on T, ro and £, should be described by scaling as: 

0(T, r , £) = e-^O^e 1 / 4 , £e 1 / 4 ) . (0.53) 

The scaling function O is such to cut-off both the singularity that appear for r Q — > oo 
at £ = £a for e > and the one at £ —> oo and e = 0. The value of the exponent a could 
be guessed if one used the potential functional as an ordinary field theory beyond the 
Kac limit for which it has been derived. Simple rescaling of the cubic action (0.43) 
leads then to a = > implying that the cross-over from Mean-Field MCT behavior 

d 

to non-critical behavior happens on length scales of the order of £\ = r 6 ~ d . The same 
conclusion could be reached analysing the behavior of the barrier for T < and 
£ < £d, which scales as r^B(£,T) ~ r$£ d ~ & {l — £/£&)■ This analysis, however, relies 
on a continuum approximation of the potential W[p(x)], which for finite ro neglects 
spatial fluctuations of the quenched disorder and heterogeneities in the reference con- 
figuration ct*- ' . It has been recently observed that finite size scaling around mean-field 
critical points in disordered systems can be non-trivial due to effective fluctuations of 
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the critical temperature (73; 74). Analogously here we can expect that effective local 
fluctuations of the critical temperatures can affect FRS. A criterion for the validity of 
naive FRS can be obtained in analogy with the analysis of FSS in disordered systems 
which in turn is an application of Harris criterion (75) for relevance of disorder in 
phase transitions. 

Let us rewrite the naive FRS of the barrier for temperatures slightly smaller than 
Td, which for I — > oo we write as 

r#Boo(T) ~ i$\ £ \"**(.<ic-<t) (o.54) 

In our case, vmf = 1/4 and d c = 6; the form (0.54) offers the possibility of studying 
more general cases. The scaling (0.54) gives a characteristic cross-over length t\ ~ 

d 

r dc ~ d beyond which deviations from mean-field theory are to be expected (remember 
that lengths are measured in units of ro). One can expect violations of this naive 
FRS when, on scale l\ the typical fluctuations of the critical temperature are larger 
then the average deviation |e|. The critical temperature fluctuations, related to local 
disorder fluctuations, can be expected to be of the order (ro(i)~ d / 2 . Applying the 

_^ d j d d c 

above criterion, we find that naive FRS is violated if r M " MF < r 2 {d - d) i.e. if 
d c VMF < 2. It is notable that this relation coincides with the criterion for the validity 
of naive FSS for disordered systems above the upper critical dimension. In our case 
d c VMF = 3/2 and so that naive scaling cannot not be expected to hold. In order to 
understand the nature of the cross-over for finite ro a deep analysis of the effect of 
critical temperature fluctuations is needed (74). This goes beyond the scope of the 
present review and won't be attempted here. 

The second point we would like to discuss is the behavior of the barrier close 
to T c . The Kac limit prediction is that the configurational entropy reduction in a 
finite volume of linear size I > l c is proportional to the surface, STi£ d = — Y£ d_1 , 
which, as observed, leads to a modified Vogel-Fulcher law for the relaxation time 
t ~ exp(rg(T — Tc) - ^ -1 )). Of course this does not necessarily imply an ideal glass 
transition if ro is finite. It just tells that for growing ro and fixed temperature T < T c 
the relaxation time grows faster then exp(ro£>) for any B > 0. The result is compatible 
both with the absence of a singularity at finite temperature or with a singularity that 
could be weaker than the one found in the Kac limit. In phenomcnological RFOT it 
is supposed that the interface exponent is renormalized by the effect of fluctuations, 
yielding ST,£ d = —Y£ e with 6 < d — 1. In rcf. (1), scaling arguments in favor of the 
value 9 = d/2 have been put forward, which leads to the canonical Vogel-Fulcher 
form in all dimensions d > 2. Correspondingly, the static length £ c would behave as 
(T — T c )~ d / 2 rather than (T — T c ) _1 as found for r — > oo. If these suggestions are 
correct, for large but finite ro there should be a cross-over between the two regimes. 
A Ginzburg criterion involving ro as well as T — T c should govern this cross-over 6 . 
Differently from the cross-over close to Td the formalism we have used does not seem 
to suggest a mechanism for this cut-off. 

6 Note that if we ask here that the fluctuations of the critical temperature on a correlated volume 
are smaller than T — T c we get the condition (T — T c ) 2> 7 — ^d/? which is always verified for 

l c ~ (T-Tc)" 1 . 
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0.5 Conclusions 

This chapter concerns the inclusion of spatial aspects in glassy theory. In particular the 
emphasis has been put on recently proposed measures of correlations -the point-to-set 
correlations- that involve an infinite number of variables. 

In the first part of the chapter we reviewed exact general bounds that relate the 
growth of a relaxation time to the growth of a correlation length. We first introduced 
with a progressive mathematical rigor the idea of point-to-set correlation functions 
and the associated correlation lengths, and explained to which extent this definition 
allows to reconcile the phenomenology of glasses with the intuitive association between 
growing time and length scales. These bounds do not spoil the interest of studying 
PS correlations in glassy systems. The lower and upper bounds have different forms, 
which correspond to different physical mechanisms of relaxation. The lower-bound is, 
up to a non-trivial exponent, of the critical dynamics type, while the upper-bound has 
the form of activated dynamics. 

In the second part of this chapter we discussed this correlation length in the context 
of the mean-field picture of glassy phenomena provided by 1RSB disordered models. 
We have thus reviewed some analytical works dedicated to growing length scales in 
simplified models of glasses and their puzzling relationship with growing time scales. 
We discussed diluted random graph models and Kac like models which add new el- 
ements to the already rich behaviour of fully-connected models of the p-spin family. 
For random graph models we found that the divergence of the point-to-set length 
accompanies the one of the relaxation time at the dynamical transition. 

The study of the Kac model reveals a rich scenario with two relevant correlation 
lengths, the mosaic length below which a system behaves thermodynamically as an 
ideal glass, and a dynamical length quantifying the typical extension of relaxation 
modes which do not require activation. At the dynamical transition the former remain 
finite, coherently with the fact that for any finite interaction range the transition is 
rounded-off, while the latter diverges. 

Let us close this discussion with a few open questions. 

The upper bound in Eq. (0.8) implies that the divergence of the equilibrium corre- 
lation time is necessarily accompanied by a divergence of the point-to-sct correlation 
length, even though with a possibly much slower form of the divergence. Both in di- 
luted random graph models, and in spin models in the Kac limit, the divergence of the 
point-to-set correlation length at Td cannot be detected by the study of static n-point 
correlation function, for any finite n. Hence in these models there is a true separation 
between the dynamic transition and the thermodynamic one, the free-energy having a 
singularity only at a lower temperature T c . An important open question is to determine 
whether this phenomenon of growth of the point-to-sct correlation function without 
a trace in the two-point function can persist in finite-dimensional models or if this is 
an artifact of mean-field models. In the latter case it would mean that bounds of the 
form (0.8) holds with £ being replaced by a standard two-point correlation length (and 
different constants C1.2). This should be important from a theoretical point of view, 
yet would not contradict the experimental and numerical situation. On the range of 
correlation time where measurements are possible the growth of the static correlations 
implied by this hypothetical extension of (0.8) could still be very weak. 
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Other questions relate to the study of finite range disordered models beyond the 
Kac limit. An accomplished theory of finite range scaling close to the dynamical tem- 
perature Td should enable to theoretically describe the crossover between mode cou- 
pling regime to activated dynamics. A related question is the connection between the 
purely static correlation length discussed here and the dynamical one extracted from 
four-points (two-time, two-location) dynamical susceptibilities. As we have discussed, 
the dynamical length found in MCT above Td can be detected in the Kac limit an- 
alyzing metastable solutions of the field equations in finite geometry. In this case, 
dynamic and static correlation length do not coincide. It would be important to better 
understand the quantitative relationship between these two characterizations of the 
correlations in a system. 
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